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Abstract 

A phenomenological framework is presented for incorporating quantum gravity motivated cor- 

O 

rections into the dynamics of spherically symmetric collapse. The effective equations are derived 

U 

from a variational principle that guarantees energy conservation and the existence of a Birkhoff 
theorem. The gravitational potential can be chosen as a function of the areal radius to yield spe- 

> 

cific non-singular static spherically symmetric solutions that generically have two horizons. For a 

(N 

specific choice of potential the effective stress energy tensor violates only the dominant energy con- 

q 

_jJ dition. The violations are maximum near the inner horizon and die off rapidly. A numerical study 

o 

of the quantum corrected collapse of a spherically symmetric scalar field in this case reveals that 

the modified gravitational potential prevents the formation of a central singularity and ultimately 

^ yields a static, mostly vacuum, spacetime with two horizons. The matter "piles up" on the inner 

H 

horizon giving rise to mass inflation at late times. The Cauchy horizon is transformed into a null, 
weak singularity, but in contrast to Einstein gravity, the absence of a central singularity renders 
this null singularity stable. 

PACS numbers: 04.25.dc, 04.70.Dy 
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INTRODUCTION 



Einstein's theory predicts its own demise in the form of cosmological and black hole 
singularities. One of the goals of quantum gravity is to provide a mathematically consistent 
mechanism for avoiding these singularities. Recently, minisuperspace models derived from 
loop quantum gravity have successfully resolved singularities in quantum cosmology pQ and 
black hole interiors [2H5]. Nonetheless, in the absence of a complete microscopic quantum 
description of black holes, it is useful to take a phenomenological (or in the more descriptive 
language of reference [B] "zoological") approach in which one tries to catalogue the possible 
non-singular semi-classical black hole spacetimes that can in principle emerge from the 
microscopic theory. This approach has a fairly long and distinguished history with much 
work focusing on the construction of non-singular static black hole spacetimes. Such metrics 
are presumed to be solutions to quantum corrected Einstein equations with an effective stress 
energy tensor on the right hand side. Bardeen [7] was one of the first to write down such 
a non-singular static black hole, but many other spacetimes have since been constructed 
[51 El E] • These solutions generically have a pair of horizons, with the inner horizon radius 
determined by the dimensionful parameter associated with the short distance behaviour. 
The conditions under which such black holes can exist were clarified in the mid-nineties by 
Borde [10] in response to the proposal in [9] for a non-singular Schwarzschild black hole 
that satisfied the weak energy condition. Inner horizons are associated with mass inflation 
[TT] and one may ask whether the resulting spacetimes are non-singular and stable under 
perturbations (i.e. tossing in more matter). 

Other methods also exist for constructing non-singular static black holes. In [T2] for 
example, the Schwarzschild metric inside a black hole was matched to a deSitter interior 
along a spacelike surface which was interpreted as a thin transition layer between the classical 
and quantum regions. Alternatively one can solve semi-classical equations of motion in a 
quantized mini-superspace model as done recently in spherically symmetric loop quantum 
gravity [2H1]. 

A potentially more fruitful approach is to examine singularity resolution in the context 
of dynamical black hole spacetimes. One of the first models was proposed many years ago 
by Frolov and Vilkovisky jT3| (see [H] for an interesting subsequent analysis). In this model 
a trapping surface forms but the inner horizon does not reach r = 0. In this case two 
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final outcomes are possible: the apparent horizons stabilize to yield a black hole with a 
static outer horizon, or the apparent horizons meet to form a closed trapping surface. In 
the latter case there is no true event horizon, although if the outer horizon is stable for 
long enough it will for all practical purposes mimic an event horizon over the relevant time 
scale. This scenario has recently been advocated by Hayward [15] . and explicitly realized 
[16] in a numerical calculation of spherically symmetric scalar field collapse using flat slice, 
or Painleve-Gullstrand (PG) coordinates (see [17] for the corresponding purely classical 
numerical analysis in flat slice coordinates). In this model, loop quantum gravity motivated 
corrections to the gravitational potential were incorporated directly into the scalar field 
equations of motion. The effect of these corrections on Choptuik scaling were analyzed, 
confirming the existence of a mass gap in the quantum corrected observed by 

Husain [18] . The choice of PG coordinates allowed the equations to be evolved beyond initial 
horizon formation thereby mapping out a large portion of the resulting spacetime, including 
the formation of trapping surfaces. For macroscopic black hole formation, two interesting 
features emerged: first the inner horizon was repelled by the repulsive core and no singularity 
formed; secondly, the modified equations of motion were non-conservative which allowed the 
outer horizon to shrink and the inner horizon to grow until they annihilated to form a closed 
trapping region. 

The purpose of the present paper is twofold. First, we present a systematic, general 
method for studying the formation of non-singular static black holes using dynamical equa- 
tions that are derived from a variational principle. A free function in the lagrangian allows 
us to specify the effective short distance behaviour of the gravitational potential which in 
turn determines the form of the vacuum solution. We look at a particular potential cho- 
sen to yield the Poisson-Israel [8] metric as a static solution. A straightforward calculation 
of the effective stress energy tensor reveals that the weak energy condition is generically 
satisfied, but the dominant energy condition (DEC) is violated outside the inner horizon. 
The violations die of asymptotically so that they are negligable near the outer horizon of 
a macroscopic black hole. For microscopic black holes these violations persist at the outer 
horizon and beyond. Since the violations of the DEC allow matter to propagate outside the 
light cone, this could contribute to the evaporation of microscopic black holes by increasing 
the rate of Hawking radiation or providing a new mechanism for mass loss. 

Second, we give the results of a numerical study of the gravitational collapse of a massless 
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scalar field using the same choice of potential as described above. The central singularity 
is indeed avoided. Despite violations of the DEC no energy escapes the outer horizon and 
there is no mass loss, thus two static horizons eventually form. The matter "piles up" at the 
inner horizon, and as expected from previous studies [H] of mass inflation, the mass function 
grows exponentially along this (nearly) null surface. The interior near the origin is flat, and 
the exterior spacetime settles down to the non-singular vacuum solution everywhere except 
along the inner horizon. 

The paper is organized as follows. In Section 2 we review several static, vacuum space- 
times that describe non-singular black holes, including the Bardeen metric and the Poisson- 
Israel metric. Section 3 describes the general action principle that we use to modify the 
gravitational potential. We write down the static non-singular vacuum solutions and derive 
the effective stress energy tensor. We also present the dynamical equations for spherically 
symmetric scalar field collapse in PG coordinates. The gravitational self-interaction of the 
scalar field is specified in terms of a single free function (the quantum corrected graviational 
potential) of the areal radius r. Section 4 reviews the numerical method for evolving the 
dynamical equations while Section 5 presents the results. We conclude with a summary and 
some speculation. 



STATIC NON-SINGULAR BLACK HOLES 



We assume that such solutions can generically be put in the form: 

ds U = - (i - -jj^y) dt * + (i - -jj^y) dr2 + M0i 2 dn a , (i) 

which is the case, for example, in the quantum corrected black hole studied in [5|. Here I is 
a parameter with dimensions of length, usually taken to be the Planck length, introduced 
for dimensional considerations. 

In general, the function j(r) will contain a parameter /i of dimension length and should 
go to r as r — > oo. Moreover, if the core is to be regular l/j(r) should vanish at least as fast 
as r 2 as r — > 0, in which case the geometry near the origin will be de-Sitter like. With these 
conditions we see that l/j(r) goes to at r = and r = oo so that there will generically be 
either no horizons, or two horizons, with a mass gap determined by the shape of j(r) and 
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the value of /i. For example, the Bardeen black hole [7] has: 

(r 2 + (j, 



On the other hand Poisson and Israel [8] used semi-classical considerations to argue for a 
non-singular static metric with: 

Their argument suggested the identification: /i 3 = a 2 M a d m , where a 2 is the quantum cut-off 
scale and M a d m is the ADM mass. This had the advantage that the size of the inner horizon 
did not go to zero as the black hole mass increased. 

In the following we will present a lagrangian formalism for deriving gravitational equations 
that produce solutions with arbitrary (but fixed) j and allows coupling to matter so that 
one can examine the dynamical formation of non-singular black holes. 



LAGRANGIAN FORMULATION AND EQUATIONS OF MOTION 



Our starting point is the action for generic dilaton gravity in two dimensional spacetime: 



1 

2G 



d 2 xy/=g (f)R(g) + S M [g, 0, ip\ + 



i 2 ) 



(4) 



where G is the 2-d gravitational constant. The dilaton <p and metric g are geometrical 
variables and Sm is a matter action that will be specified below. The most general solution 
can be written in the form: 



ds 2 



i 



21GM\ 



dt 2 



w ) 



(5) 



j{4>) ) \ ./( 

where M. = M.(t,r) is a generalized Misner-Sharpe mass function that is related to the 
traditional definition [20] by A^trad — r A4/j(0). The vacuum theory with a constant mass 
M satisfies a Birkhoff theorem [2TJ . j(<p) is related to the dilaton potential by: 

dj 



d(fi 



v{4>) 



and r = r(0) is determined by: 



dr = I 



dcj) 



(6) 



(7) 
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(n-l)\l 



n 



This action describes a wide class of theories containing black hole solutions, including 
the spherically symmetric sector of D-dimensional Einstein gravity. The latter requires the 
identifications: 

167rG (n+2) n 

2G = 8(n-l)z»K (8) 

n fr\ n , , 

(9) 

l/n 

y W = ^[^)) (10) 

where G( n+2 ) is the D-dimensional Newton's constant, r is the radius of a rotational invariant 
two-sphere, and 

"-m*® (12) 

is the volume of the n-dimensional unit sphere. The physical D-dimensional metric is: 

1 



dS U = + r2 (^) dn (n) > ( 13 ) 



where dQ 2 n ^ is the volume element of the unit n-sphere. We have chosen to use h(r) as for 
D-dimensional Einstein gravity, though one is free to choose a different form. It is straight- 
forward to verify that with the above identifications, the following D = n + 2 dimensional 
metric is the correspondng Tangerlini-Schwarzschild solution: 



2l D ~ 2 GM\ / 2l D ~ 2 GM 



ds 2 (D) = -[l- rD _ 3 )dt*+[l- rD _ 3 j dr 2 + r 2 <if!f„, . (14) 

In the following we will focus on n = 2, i.e. 4-dimensional spherically symmetric gravity. 
The generalization to higher dimensions will be straightforward. 

We propose that the action Q be used to generate the dynamical equations of motion 
for spherically symmetric collapse in 4-dimensions with quantum corrected gravitational 
potential. In particular, one defines the physical 4-d metric in terms of 2-d quantities by 



(13), first choosing the function j(r) and then specifying <p(r) and V(0) in the action to give 



this j(r). That is, one requires: 



ld(j) = j(r)dr ; (15) 

VU) = d ^ = dT = 1 ^ (16) 
d(f) dr d(f) j(r) dr 



For example, for the Poisson-Israel metric ([3]), which will be the focus of our numerical 
studies, one obtains: 

m = f dr^±p- = % - 4- ; (17) 
rv ; / 2lr 2 41 2lr ' V ; 



= "w ■ < 18 > 



l{r 3 - 2/i 3 ) 
r(r 3 + fi 

The above procedure yields a vacuum solution that corresponds to the static spherically 
symmetric 4-d metrics of precisely the form [I]). Note that both and V(<p) go to minus 
infinity as r — > 0. This behaviour is generic for any function with asymptotic form j(r) ~ r~ a 
with a > 1. Only j(r) enters the equations of motion in the subsequent analysis, so these 
divergences do not affect the resulting dynamics. 

Once the gravitational part of the action has been chosen, one can study the quantum 
corrected dynamics of spherically symmetric collapse by specifying the matter lagrangian. 
For simplicity we consider a massless scalar field, but there is no obvious impediment to 
study other forms of matter, including electromagnetic and Yang-Mills fields. Our matter 
lagrangian is: 

L M = J d 2 xr\x)^-g\V^\ 2 , (19) 

where we have made explicit the fact that r(x) is a scalar function of the coordinates, V 
refers to a 2-d gradient, and we have been consistent in using h(r) = r 2 in this lagrangian 
as well. 



Effective Stress Energy Tensor 



Given the 4-d metric (13) it is straightforward to calculate the 4-d Einstein tensor and 



write the equations of motion in 4-d form. The result is: 

G("„> = 8fG< 4 >^ (v„VV„V " lC IWf) + It ; (20) 

G S' = -y G<4) 2/7M |w|2+7 ^' (21) 

In the above [i, v — 0, 1 refer to the spherically symmetric coordinates (t, x), while i, j = 2, 3 
refer to the angular coordinates. The modified action gives rise to two effects. First, there 
is an effective r-dependent gravitational constant: 

G eff = G (4) _J_ = ^(4) (22) 

2lj(r) r 3 + /i 3 ' v ; 



which vanishes as r — > 0. Secondly, there is an effective stress energy tensor of the form: 

_ 2G^M M _ 1\ (4) (4) 

2^M . - , . (2 1) 




where the ' denotes differentiation with respect to r. This effective stress tensor is only 
non-vanishing when the mass function is non-zero. 

The mass function determines the location of apparent horizons via: 

1 - — — = . 25 

j(r) 

If one uses the traditional definition of the Misner-Sharpe mass function then the horizon 
location is given by: 

1 - 2G * M «>* = o . (26) 
r 

Since r/j(r) is a monotonic function for the class of theories we wish to consider, as long 
as the matter action is standard, violations of the energy conditions can only come from the 
effective stress tensor. Its simple form allows a detailed analysis. In fact, one can verify the 
following inequalities are required by the corresponding energy conditions: 

p < Weak Energy Condition (WEC) ; 
(3 - p > Null (NEC) and Strong (SEC) Energy Conditions ; 
\p\ - \(3\ > Dominant Energy Condition (DEC) . (27) 



As shown in Figs. [jj[3j the WEC, NEC and SEC are satisfied for the j(r) that we are 
using, but the DEC is violated outside the inner horizon. Notice that l/j{r) has a positive 
slope for r < 2 1//3 /i, which represents the extent of the repulsive core. Correspondingly, the 
DEC violations begin at this point and drop off rapidly as r -6 , so that they are not seen 
near the exterior of large black holes. However for microscopic black holes the violations 
can extend beyond the outer apparent horizon. Since the DEC dictates that mass-energy 
cannot flow faster than light, these violations could play a role in describing the evaporation 
of microscopic black holes within the theory. 



Solving the horizon condition (25) we find that the inner and outer horizons coincide at 
r = 2 l l z jjb giving the size of the smallest black hole allowed to form. This signifies a mass 
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FIG. 1: A plot of the weak energy condition for \i = 1 and M = 5. 




FIG. 2: A plot of the null and strong energy conditions for \i = 1 and M = 5. 



gap in the theory with a threshold mass of: 



M, 



th 



3(2) 1 / 3 /i 
4G( 4 ) 



(28) 



Hamiltonian Equations 



In order to look at spherically symmetric collapse and horizon formation it is convenient 
to use flat slice, or PG coordinates. Choosing the space coordinate x — r, the generalized 



9 




FIG. 3: A plot of the dominant energy conditon for \i = 1 and M = 5. With these parameters, the 
inner horizon is located at ~ 0.225 while the outer horizon is at r h ~ 10. The curve crosses 
the axis at r = 2 1 / 3 . 



PG metric is given by 



ds 2 ^ 



-a 2 dt 2 + [ dr + \ l ^^— adt ) + r 2 dtt 2 , 



(29) 



where a is the lapse function which scales time across a spatial slice. The gauge fixing 
procedure that determines the reduced Hamiltonian and equations of motion for the scalar 
field were derived in [TTJ . One first makes the gauge choice 



which yields the consistency condition on the lapse and shift functions: 



aGIi p = Ncj)'. 



(30) 



(31) 



The equations can be put in a more transparent form using the canonical transformation 



X := e p 



P := e- p GIi p , 



(32) 
(33) 



where U p is the momentum conjugate to p. P is then conjugate to X with their Poisson 
bracket being 

{X(x),P(y)} = G5(x,y). (34) 
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The resulting Hamiltonian is 



dr a I --^-rM + Qm + I 



+ ldr[ W) M 



(35) 



where 



M 1 I P- 

2G 



1 / % 

2 U(0) 



X 2 I 2 



(36) 
(37) 



Ai is the Misner-Sharpe mass and approaches a constant at spatial infinity where it is equal 
to the ADM mass of the solution. Qm is the energy density of the scalar field. 
To completely fix the gauge we choose: 



(38) 



This condition, along with Eqs.(36) and (30), implies that P 2 = 2GM./1. These gauge 
conditions produce the non-static generalization of PG coordinates, as can be seen by using 



(31) and the gauge conditions in the case of a vacuum to reduce the line element to PG 



form [T7] : 



ds 2 =j(<P) 



,, , , 2GMI , 
-at + | dr + \ I — r-^—dt 



(39) 



More importantly, the spatial slices are regular across apparent horizons that form during 
the evolution. 

We are now able to write the equations of motion for the scalar field in fully reduced 
form: 



. ly/2GM/W , 

W = & I 7=^ 1" 



K4>) ) ' 



n, 



(40) 
(41) 



where a and M. are the solutions to 



M 1 = Q M + lip% 



'2GMI 

1W 



(42) 
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a H =cr = 0. 43 



The right hand side of (42) defines the instantaneous mass density of the configuration. 
It has the expected contribution from the energy density of the scalar and an additional 
coupling term that corresponds to the contribution from its self gravity. 

The above equations need to be supplemented by boundary conditions for M. and a. 
Without loss of generality we choose a = 1 at r = 0. A change in this value corresponds to 
a trivial rescaling of the time coordinate. For future reference we note that this condition 
implies that PG time t and proper time s are equivalent at the origin. We also fix M. = at 
r = 0, which guarantees that the metric is flat in the neighbourhood of the origin. One way 
to think about this is that with this identification we have eliminated the single gravitational 
degree of freedom associated with the spherically symmetric gravitational field [19J so that 
the gravitational mass is provided totally by the mass content of the spacetime. 

NUMERICAL METHODS 

The first step in the iteration process is to specify initial if) and 11^ configurations. We 
work with two forms: 

2" 



if) = At exp 



r — tq 
B 



(44) 



rj, = Atanh(^-^j, (45) 

where A, B and r are the parameters which dictate the initial mass function. Note that 
since the mass density depends on if)', the tanh form has one mass peak while the gaussian 
data has two. For both cases we define an initial standing wave by choosing EU^(r, 0) = 0. 
The results were identical for both forms of initial data. 

In order to have a fine grid spacing where the dynamics require a more accurate inte- 
gration, we use a variable r-spacing that remains unchanged throughout the evolution. A 
typical r-lattice would use a spacing of Ar(r) = 10^ 3 from the origin until some point be- 
yond the inner horizon, at which point the lattice increased by one percent each step until 
reaching a maximum of 10 -2 . With this maximum step, the lattice would be extended as 
far as necessary to contain all of the mass within the system (numerical errors occur when 
mass leaves the grid). 
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We use an adaptive time step At(t) refinement according to the minimum found across 
the spacial slice with the condition 

At(t) = MIN r | ^Ar(r) J , (46) 

where ^ is the inverse of the local speed of an ingoing null geodesic. This provides stability 
by preventing information from moving over too many r-points in a single time step. 

For derivatives, we experimented with finite differences and cubic splines finding the same 
qualitative behaviour using both, which gives some indication of reliability. In the end we 
chose to use finite differences exclusively since they require less computational time. The 
variable lattice was problematic for fourth order finite differences, so we used second order 
differences. 

After specifying intial if) and ILy,, fourth order Runge-Kutta (RK4) methods are used to 



find Ai and a by integrating (42) and (43) across the spatial slice. These values are then 
entered into the right hand sides of (40 ) and (41 ) to find ip and at the next time step with 
RK4. This process is repeated until the matter disperses to leave behind a flat spacetime, 
or in the case of black hole formation, for as long as the code remained stable. In the latter 
case, the mass density was seen to bounce off the origin and move outward toward the inner 
horizon. As the mass moves closer and closer to the inner horizon, instability eventually 
sets in. However, the code always runs for long enough to extract the relevant information. 

Convergence of our code is apparent since the results were not affected by increasing the 
resolution beyond 10 -2 as long as the black hole size on formation was sufficiently large 
enough compared to the grid size. To ensure the stability of our code, we monitor the ADM 
mass for each run and find it to remain constant within ~ 0.1%. 



RESULTS 



With the quantum scale parameter /i set to zero we find results identical to those in 
|17j . If the mass profile on the initial time slice evolves such that Ai(r = 2 1/,3 /i) < My, at 
all times, then no horizons form and the mass is free to disperse to infinity leaving behind 
a flat spacetime. If the threshold mass is reached, black hole formation begins with the 
appearance of a single horizon that splits into two as the evolution continues. The outer 
horizon continues to grow as mass falls through it while the inner horizon moves toward the 
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origin. The horizon paths define a two dimensional trapping surface in spacetime within 
which all null geodesies point toward the origin, illustrating several features expected in 
theory [T5]. Specifically the trapping horizon allows for multiple apparent horizons during 
intermediate phases of the collapse which annihilate in pairs. It necessarily terminates at 
the singularity at one end, and at null infinity (the horizon) at the other. 

For non-zero jx we find the mass of the smallest black hole to be in good 



agreement with the analytic calculation for static black holes (28). The outer horizon dy- 



namics are similar to the classical case, however the repulsive core causes the inner horizon 
to slow as it moves toward r ~ fi, eventually settling down to a constant value after all of the 
mass is inside the inner horizon. The mass continues moving inward toward the origin, but 
is now free to move through the origin and back out toward the inner horizon. As this mass 
nears the inner horizon, it continually slows so that it does not enter the trapping region in 
finite proper time. See Fig. [4] for a spacetime diagram showing the horizon locations and 
null geodesies throughout the evolution 




space 



FIG. 4: A space-time diagram showing the trapping surface (red, dashed line) and null geodesic 
(black, solid lines). Notice all "outward" moving null geodesies inside the black hole converge at 
the inner horizon. 



To investigate mass inflation, we solve for outgoing null geodesies by integrating: 

-i 
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We then define a coordinate v along each outgoing null geodesic by: 

-l 

hmu\ 

dv = adt + 1 + 



J 



j dr . 



(48) 



Note that since we are not including the appropriate integrating factor the above does not 
define the coordinate transformation to Vaidya time. It does, however, define an affine 
parameter along each outgoing null geodesic and is sufficient to examine the issue of mass 
inflation. In particular, previous studies of classical black holes possessing an inner horizon 
(i.e. charged or rotating black holes) have found that along outgoing null geodesies origi- 
nating near the origin, the mass function increased exponentially with Vaidya time, forming 
a weak [25J mass inflation singularity near the Cauchy horizon [Til [221 [23]. In these sce- 
narios the Cauchy horizon eventually contracts to the origin and forms a strong, spacelike 
singularity. 

The plots of ln(.M) vs. v in Fig. [5] confirm mass inflation within the present quantum 
model. However in this case the inner horizon settles down to a constant radius, thereby 
removing the central strong singularity. A drawback of this model is that the final location 
of the inner horizon tends to zero for fixed fi with increasing mass, inevitabley forcing the 
mass toward the origin for macroscopic black holes. 




FIG. 5: Plot of ln(M) vs. v for a series of null geodesies leaving the origin after black hole 
formation. 



Note that the dynamical equations violate the DEC similarily to the vacuum solution, 
however this was not seen to allow for any black hole mass loss. 
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CONCLUSION 



We have developed a Hamiltonian formulation of spherically symmetric gravitational col- 
lapse in the presence of quantum corrections. The form of the correction is freely specifiable 
and provides a computational laboratory for observing gravitational collapse in the presence 
of effective potentials arising from modifications to general relativity. 

With the choice of j(r) studied herein, evolution of the dynamical equations results in 
either a flat spacetime or a black hole possessing two horizons, with the mass piling up 
along the inner horizon and forming a weak, null singularity. This family of corrections also 
permits a stable, vacuum solution. The collapse scenario and the static solution possess an 
equivalent mass gap that is proportional to the quantum scale parameter. 

The modified potential was found to violate the DEC, which can allow for mass energy 
to flow faster than light. However, this effect was not observed in the numerical evolution; 
mass that entered the black hole was found to be trapped inside the inner horizon. A mass 
inflation singularity was found to form along the inner horizon. This horizon was stationary 
at late times rather than contracting to the origin as in classical Kerr or Reissner-Nordstrom 
black holes. 

An improvement to this model would be to consider some form of j(r) that scales the 
quantum parameter fi so that the inner horizon does not go to zero for macroscopic black 
holes, similarly to the static solution studied by Israel and Poisson. Also, it would be quite 
interesting to incorporate a radiation term in the Hamiltonian to allow for evaporation which 
would likely lead to a closed trapping region as in [15l [16] . In this scenario one wonders how 
mass inflation would play out as the black hole shrinks and eventually vanishes. 
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